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Enhancing the polarization of nuclear spins surrounding the nitrogen-vacancy (NV) center in diamond has 
attracted widespread attention recently due to its various applications. Here we present an analytical theory and 
comprehensive understanding on how to optimize the dynamic nuclear polarization by an optically pumped NV 
center near the ground state level anticrossing. Our results not only provide a parameter-free explanation and 
a clearly physics picture for the recently observed polarization dependence on the magnetic field for strongly 
coupled *^C nuclei [H. J. Wang et al., Nat. Commun. 4, 1 (2013)], but also demonstrate the possibility to 
strongly polarize weakly coupled nuclei under weak optical pumping and suitably chosen magnetic field. 

This allows sensitive magnetic control of the nuclear spin polarization for NMR applications and significant 
suppression of the '^C nuclear spin noise to prolong the NV spin coherence time. 


INTRODUCTION 

The atomic nuclear spins are central elements for NMR 
and magnetic resonance imaging! and promising candidates 
for storing and manipulating long-lived quantum information! 
due to their long coherence time. However, the tiny magnetic 
moment of the nuclear spins makes them completely random 
in thermal equilibrium, even in a strong magnetic field and 
at low temperature. This poses severe limitations on their 
applications. The dynamic nuclear polarization (DNP) tech¬ 
nique can bypass this limitation by transferring the electron 
spin polarization to the nuclear spins via the hyperfine inter¬ 
action (HFl), but efficient DNP is usually prohibited at room 
temperature. 

An exception is the negatively charged nitrogen-vacancy 
(NV) center— in diamond, which has an optically polariz¬ 
able spin-1 electronic ground state with a long coherence 
time!, allowing DNP at room temperature!*!. This prospect 
has attracted widespread interest due to its potential applica¬ 
tions in room-temperature NMR, magnetic resonance imag¬ 
ing and magnetometryi!, electron-nuclear hybrid quantum 
register——, and electron spin coherence protection by sup¬ 
pressing the nuclear spin noise!!. In addition to the remark¬ 
able success in coherently driving spectrally resolved transi¬ 
tions to initialize, manipulate, and readout up to three strongly 
coupled nuclear spins!!*!!^— , there are intense activities aim¬ 
ing to enhance the polarization of many nuclear spins via 
dissipative spin transfer from the NV to the nuclear spins. 
To overcome the large energy mismatch for resonant spin 
transfer, various strategies have been explored, e.g., tuning 
the NV spin near the excited state level anticrossing!*!!^— 
or ground state level anticrossing (GSLAC)!*!!*!!, driving the 
NV-nuclear spins into Hartman-Hahn resonance!!*!! or selec¬ 
tively driving certain spectrally resolved transitions between 
hyperfine-mixed states under optical illumination!!*!!. Suc¬ 
cessful polarization of bulk nuclear spins in diamond have 


dramatically enhanced the NMR signal by up to five orders 
of magnitudes!!*!! and significantly prolonged the NV spin 
coherence time!!*!!. 

In particular, near NV excited state level anticrossing, al¬ 
most complete polarization has been achieved for '^N (or '^N) 
and nuclei in the first shell of the vacancy!*!!^i!!*!!. Re¬ 
cently, Wang et al?^ exploited the GSLAC to achieve near 
complete polarization of first-shell nuclei and revealed 
multiple polarization sign reversals over a narrow range (a 
few mT) of magnetic field. This behavior has been attributed 
to the anisotropic HFI and could allow sensitive magnetic 
control of nuclear spin polarization, but a clear under¬ 
standing remains absent. Furthermore, in most of the existing 
works, only a few strongly coupled nuclear spins (HFl » 200 
kHz) are significantly polarized via direct spin transfer from 
the NV center, while many weakly coupled nuclear spins are 
only slightly polarized via nuclear spin diffusion. Enhancing 
the polarization of these weakly coupled nuclear spins could 
further improve NMR and magnetic resonance imaging!!*!!*!! 
and prolong the NV spin coherence time!!*!!. 

In this paper we present an analytical formula and a com¬ 
prehensive understanding on how to optimize the DNP by an 
optically pumped NV center near the GSLAC at ambient tem¬ 
perature. Our results provide a parameter-free explanation and 
a clear physics picture for the experimentally observed mag¬ 
netic field dependence of nuclear polarization!!. More 
importantly, we demonstrate the possibility to greatly enhance 
this magnetic field dependence and strongly polarize weakly 
coupled nuclei (with HFI down to ~ 1 kHz) via direct, 
resonant spin transfer under suitable conditions. First, we 
introduce our model and an intuitive picture for manipulat¬ 
ing the DNP with an optically pumped NV center. Then we 
present our theory and analytical formula for the DNP of a 
single nuclear spin, which reveals the possibility of strongly 
polarizing weakly coupled nuclear spins under weak optical 
pumping and fine-tuned magnetic field. Finally, we demon¬ 
strate this possibility for multiple nuclear spins by providing 
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FIG. 1. (a) NV states at ambient temperature responsible for DNP 
under optical pumping, (b) Reduced five-level model for DNP near 
the ground state anticrossing. 


a multi-spin DNP theory and its numerical solution. 


RESULTS 

Model and intuitive physics picture 

First we briefly introduce our theoretical model, consisting 
of a negatively charged NV center coupled to many surround¬ 
ing nuclear spins at ambient temperature. The NV center has 
a ground state triplet 1+1^) and |0g) separated by zero-field 
splitting Dgs = 2.87 GHz and an excited state triplet | + Ig) 
and |0g) separated by zero-field splitting = 1.41 GHz22. 
In a magnetic held B along the N-V axis (z axis), the elec¬ 
tron Zeeman splitting yeB with = 28.025 GHz/T can¬ 
cels the ground state zero-field splitting at a critical magnetic 
held Dgslje ~ 102 mT, leading to GSLAC between |0g) and 
I - Ig). The GSLAC reduces the energy mismatch for the 
NV-nuclei flip-flop dynamics and enables NV-induced DNP 
through their HFI i^HF = where F, s Sg-Ag /H-Sg-Ag ,■ 

is the electron Knight held coupled to the /th nucleus, and 
Sg (Sg) is the NV ground (excited) state spin. The HFI ten¬ 
sors Agj and A^j are well established for the on-site nitrogen 
nucleus^ and the nearest nuclear spins in the first shell 
of the vacanc y For other nuclei, especially those 

> 6.3 A away from the NV center, very little information is 
available ^^i^^~ — and dipolar HFI is usually assume d^^'^^'^^^ . 

Next we provide an intuitive physics picture for engineer¬ 
ing the DNP by controlling the NV center near the GSLAC. 
For brevity we focus on one nuclear spin-/ (e.g., 7=1/2 for 
the '^C or '^N nucleus and 7=1 for the '^N nucleus) and 
drop the nuclear spin index i. To describe the NV-nucleus 
flip-flop, we decompose the HFI into the longitudinal part 
Fj/j. that conserves the nuclear spin and the transverse part 
(F+I- + h.c.)f2 that induces NV-nuclear spin flip-flop dy¬ 
namics. As we always work near the GSLAC, the nuclear 
spin flip induced by the off-resonant excited state HFI is very 
small (to be discussed shortly). Since | - Ig) is nearly de¬ 
generate with the NV steady state |0g) under optical pumping, 
the NV-nucleus flip-flop is dominated by the pro¬ 


cess |0g) ® \m) —> I - Ig) ® |m -F 1) and the process 

|0g) ® \m + 1) —> I - Ig) ® \m), where {\m)} are nuclear spin 
Zeeman eigenstates. The corresponding energy mismatches 
of these two processes are 

= A -F y^B (m -F l)Ag^^2, 

Am(^m+\ = A yj\fB ITlAgg^^, 

respectively, where the first term A s Dgs - yeB is the | - Ig)- 
|0g) energy separation, the second term ±yNB comes from the 
nuclear spin Zeeman term, while the last term comes from the 
longitudinal HFI fJ^ ^ -(Ag + Ae^zz^-ie-ie)h- The 

energy mismatch difference 

An = — 2yNB Ag^^^ ( 1 ) 

provide the first essential ingredient for engineering the DNP 
by selectively driving one process into resonance. 

The second essential ingredient for engineering the dissi¬ 
pative DNP is optical pumping, which connects the NV cen¬ 
ter to the strongly dissipative bath of vacuum electromagnetic 
fluctuation and makes the NV center itself a dissipative, non¬ 
equilibrium bath: it quickly establishes the unique steady state 
|0g) whatever the initial state is (commonly known as opti¬ 
cal initialization). When optical pumping is not too weak, 
such that the NV dissipation (more precisely the NV opti¬ 
cal initialization) is much faster than nuclear spin dissipation 
(more precisely the DNP process), the NV center can be re¬ 
garded as always in its steady state |0g), i.e., it becomes a 
Markovian bath and induces DNP by driving the two pro¬ 
cesses The resonance linewidth of both processes 

are determined by the optically induced NV ground state level 
broadening. Therefore, under sufficiently weak (but not too 
weak) optical pumping such that the linewidth of each reso¬ 
nance is much narrower than the distance lA^fl between the 
two resonances, we can fine tune the magnetic held to drive 
into resonance, while keeping 
(o‘-ij,0g7+) off-resonance to achieve strong positive (negative) 
nuclear polarization. When the magnetic field is swept from 
resonance to resonance, the nuclear po¬ 

larization changes from positive values to negative values over 
a magnetic field range ~ \A.N\lye- For first-shell nuclei, 
IAatI is dominated by |Ag ;,j| » 130 MHz. This gives a simple 
explanation to the experimentally observed polarization sign 
reversal over a few mT^'*. For weakly coupled nuclei, 
IAatI ~ 2\yNB\ X 2.2 MHz is greatly reduced, corresponding to 
a much more sensitive dependence of the polarization on 
the magnetic field over a range ~ 0.1 mT. Below we demon¬ 
strate this intuitive idea and provide a complete picture for 
the DNP of a single nuclear spin-7 before going into the more 
complicated case of multi-nuclei DNP. 


DNP theory of single nuclear spin 

Under optical pumping from the NV ground orbital |g) 
to the excited orbital \e), seven NV states are relevant, in¬ 
cluding the ground state triplet, excited state triplet, and a 
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metastable singlet state \S) [see Fig. [TJa)]- The NV Hamil¬ 
tonian consists of the orbital part a>Q&e,e, the optical pumping 
term Hdt) - (Q.Rl2)e~'‘^'&e,g + h.c., the ground state triplet 
part Hgs s Dgs5g j + yeBS and the excited state triplet part 
Hes s -I- jeBSe,z, where &i j = |/){y| is the transition 

operator. In the rotating frame of the optical pumping, the 
NV-nucleus coupled system obeys 

p(f) = -Cnvp (0 - i{HN + F ■ I,p(f)]. ( 2 ) 

where = JnBIz is the nuclear spin Zeeman Hamilto¬ 
nian and -Cnv(-) = -*[^nv,(-)] + Yja^\Ba\(-) is the evo¬ 
lution superoperator of the NV center in the absence of the 
nuclear spin, with Hnv the NV Hamiltonian in the rotating 
frame and the last term accounting for various NV dissipa¬ 
tion channels [shown in Fig. [Ua)] in the Lindblad form 
D[L]p = {LpD - {DLp]i2), as well as the excited orbital 
pure dephasing at a rate F^ and the ground state spin dephas¬ 
ing at a rate jg, - 1 / 7 ’ 2 ,nv that models the hnite ground state 
coherence time r 2 ,Nv- We use the experimentally measured 
dissipation rates at room temperature; Fg = 10 ^GH 2 S^, 7 = 13 
MHz22, 71 13.3 MHziS, and 7 , = 0.56 MHzil^. Here we 

have neglected the very small leakage via intersystem cross¬ 
ing from lOe) to I + Ig), consistent with the experimentally 
reportedi^i^ high optical initialization probability ~ 96% into 
|0g). Indeed, we have verified that the nuclear polarization is 
not so sensitive to the optical initialization probability, e.g., 
upon including appreciable leakage such that the optical ini¬ 
tialization probability drops to 80%, the steady-state nu¬ 
clear polarization drops by ~ 20 %. 

The DNP can be directly obtained by solving Eq. (I2]l nu¬ 
merically. However, this approach does not provide a clear 
physics picture for the underlying DNP mechanism and is not 
suitable for searching for optimal experimental parameters to 
maximize the DNP effect. Further, for the DNP of multiple 
nuclear spins (to be discussed shortly), this approach quickly 
becomes infeasible with increasing number N of the nuclei, 
because the dimension of the Liouville space grows exponen¬ 
tially as dim(p) » {21 + 1)^^M^, where M = 7 is the number 
of relevant energy levels of the NV center, e.g., the DNP of 
N - 5 nuclear spin-l/2’s already involves dim(p) a; 50000 
and hence is computationally intensive. 

Our work is based on a recently developed microscopic 
theory^^^, applicable as long as the NV dissipation (opti¬ 
cal initialization in the present case) is much faster than the 
nuclear spin dissipation (DNP in the present case). The popu¬ 
lations lp,„} of a single nuclear spin-/ on its Zeeman sublevels 
|m) (m - -I, -I + 1, ■ ■ ■ , /) obey the rate equations 

pm ~ 'y ] ^m<—nPn ( ^ (3) 

n=m±l n=m±l 

Up to second order of the tranvserse HFI (longitudinal HFI 
treated exactly), the transition rate from \m) to \m + 1 ) is^ 

^+oo 

= ^Re dt (4) 

^ Jo 

where s (OT|/:p/±|m), the Liouville superoperator i]„ ,„(») s 
■£nv(*) - i{nFzi*‘) ~ i*‘)niFz\, and is the NV steady state in 


the rotating frame conditioned on the nuclear spin state be¬ 
ing \m), i.e., - 0 and TrNvPm = 1- The nuclear 

spin depolarization by other mechanisms such as spin-lattice 
relaxation can be described by a phenomenological depolar¬ 
ization rate 7 dep and incorporated by replacing with 

m + ^ftj7dep- 

Equation (|4| can be understood as a generalized fluctuation- 
dissipation relation^. Since the nuclear spin transition from 
\m) to \m -I-1) (|m - 1)) is induced by the HEI term F_/+ {F+IJ) 
and involves the transfer of energy jf^iB from the NV cen¬ 
ter (nuclear spin) to the nuclear spin (NV center), the corre¬ 
sponding transition rate is proportional to the fluctuation of 
the Knight field (/"+) at frequency —jnB i+y^B) in the 
NV steady state P„,. Without the longitudinal HEI we 
have Xn.m = -Cnv, so P,„ = P is the steady NV state in the 
absence of the nuclei, and Eq. (|4| reduces to the conventional 
Born-Markovian approximation. As discussed above and fur¬ 
ther elaborated below, the longitudinal HFI F^Iz could sig¬ 
nificantly modify the energy cost of the nuclear spin flip and 
hence the steady-state nuclear polarization. 

In general, the nuclear spin transition rates in Eq. 

(|4]i can be evaluated numerically (see Methods). However, the 
physical picture and global optimization of the experimental 
parameters will be greatly facilitated by explicit analytical ex¬ 
pressions, which we present below. 


Explicit analytical expressions 


Since we always consider B ~ 102 mT near |0g)-| - Ig) 
GSLAC, the dominant NV-nucleus flip-flop involves the near¬ 
resonant process |0g) |-lg)- This allows us to neglect I-Hi g) 
and I -H Ig), describe the NV center by a hve-level model [Eig. 
(lib)], and derive explicit analytical expressions for the nu¬ 
clear spin transition rates Wm±\i~m (see Methods). Eor general 
parameters, the expressions are complicated due to the pres¬ 
ence of various quantum coherence effects associated with the 
optical pumping induced | 0 g)-| 0 <;) and | - lg)-| - U) oscilla¬ 
tion. However, if we focus on weak optical pumping far from 
saturation, i.e., optical pumping rate R = ^ linewidth 

7 -H 71/2 » 26.3 MHz of |0e) I-Ig) transition, then the quan¬ 
tum coherence effects are negligible and Wm±\^,n assume the 
appealing form of a Eermi golden rule 




g.+r 


2V2 


5T)(A„ 


(5) 


where Pg - {R + y)l(2R + 7 ) is the steady-state popula¬ 
tions on |0g), b*^*(A:) s {ylji)l{x^ + 7 ^) is the Lorentzian 
shape function, F s 7^ -h /? is the resonance linewidth, and 
s e+ ■ Ag ■ CiF (e± s ± ie.y) quantifies the coeffi¬ 
cient of the term in the transverse HFI. A key fea¬ 

ture of Eq. (|5]l is that the resonance linewidth of both nuclear 
spin-flip processes is the sum of | 0 g)-| - Ig) decoherence rate 
7ip - 1 /T 2 ,nv (~ 1 kHz <K typically R) and optical pump¬ 
ing induced NV level broadening (which equals the optical 
pumping rate R). This result provides a clear understanding 
of the laser induced NV ground state level broadening previ- 
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ously observed in the cross-relaxation between the NV cen¬ 
ter and nearby nitrogen impurities^^. Further, it indicates 
that under sufficiently weak pumping, the opposite processes 
Wm+i^m and Wm^m+i Can be selectively driven into resonance 
to achieve strong nuclear spin polarization. Equation (|5]l also 
reveals that the DNP depends strongly on the form of the HFI 
tensor Aj, e.g., for the or nucleus with isotropic trans¬ 
verse HFI Ag^xx - Agjy and Ag^^y - 0, we have ++ = 0 and 
Ag_+_ = '^Ag^xx and hence complete positive polarization. 

Equation Q can be applied to a general nuclear spin-/, but 
hereafter we focus on a nuclear spin-1/2 and define W+ = 
VF-c^X and W- = VFx^|, given by Eq. 0 as 


VF± = Pg2n 


2 V2 


2 

+ Aa,/2). 


( 6 ) 


The evolution of the nuclear polarization p = 2{Ig} follows 
from p{t) = -W[p{t) - pss] as p{t) = p^^ + [p(0) - Pss]e^"'', 
so the DNP of a '^C nuclear spin is completely determined by 
two key quantities: the steady-state polarization p^^ = (VF+ - 
VF_)/(VF+ + 1T_) and the rate VF s Wy + IF_ of DNP. Up to now, 
we have neglected other nuclear spin relaxation mechanisms, 
such as spin-lattice relaxation and transverse HFI with NV 
excited states. The former occurs on the time scale ranging 
from a few seconds to tens of minute a'^i^^ . The latter occurs 
on the time scale 1/Wes, where Wes can be estimated from 
the Fermi golden rule as 2n{Ael2)^6^^^^^l^'*{D^P)Pe, e.g., for 
a '^C located at |R| > 3 A coupled to the NV excited states 
via dipolar HFI, we have 1/Wes > 1 s for /? = 0.2 MHz. 
These effects are described by nuclear depolarization W± —> 
W± -H 7depj which increases W by 27dep and decreases pss by 
a factor 1 -H 2ydep/IF, i.e., nuclear depolarization is negligible 
when W » 2ydep- 

Finally we emphasize that our rate equation theory and 
hence Eq. (|6]l are accurate when the DNP time 1 /W from 
Eq. ® » NV optical initialization time Tc ~ {{RyiKjx + 
T)]~* + T/* ■ When the optical pumping is so weak that the cal¬ 
culated DNP time from Eq. ® drops below t^, the NV center 
becomes a non-Markovian bath and the true DNP time would 
be lower bounded by ~ Tc, as we discuss below. 


DNP of strongly coupled nucleus 

To begin with, we consider the DNP of a strongly cou¬ 
pled nucleus in the first shell of the NV center under 
optical pumping near the |0g)-| - Ig) GSLAC. This configu¬ 
ration has been studied experimentally in ensembles of NV 
centers^l, which reveals an interesting polarization sign rever¬ 
sal over a narrow range of magnetic fields. This is attributed 
to the anisotropic HFI and could allow magnetic control of 
the nuclear polarization, but a quantitative understand¬ 
ing remains absent. Here we provide a clear physics picture 
and parameter-free explanation for this sign reversal using 
the experimentally measured ground state HFI tenso r'^i^^i^'i^^ 
Ag,xx = 198.6 MHz, Ag,j,^ = 123.0 MHz, Ag^,, = 129.0 MHz, 
Ag xz = Ag.zx - -21.5 MHz, and excited state HFI tensor— 
= 103.2 MHz, A,^yy = 56.7 MHz, A.^^z = 79.5 MHz, 



FIG. 2. (a) Pss of a first-shell *^C nucleus: exact solution to Eq. 

S with (orange solid line) or without (black solid line) excited 
state HFI, compared with the experimental data^i (filled squares and 
empty circles) from two different analysis methods, (b) Real-time 
evolution of nuclear polarization p(t) on VF+ resonance (A = - Aiv/2) 
or far from the resonance (A = -0.44 GHz): exact solution to Eq. (|2j 
vs. our rate equation theory [Eqs. and ©]. 


Ae,xz - Ae,zx - “32.6 MHz. To focus on the intrinsic behav¬ 
ior of DNP induced by the ground state HFI, we set ydep = 0. 
We have verified that due to the strong HFI induced NV-^^C 
mixing, the magnetic field dependence of the nuclear po¬ 
larization depends weakly on the optical pumping up to /? < 4 
MHz. In our numerical calculation, we take P - 0.4 MHz. 

First we discuss the exact solution to the Lindblad master 
equation Eq. I® without the excited state HFI. In Fig. |2ja), 
the calculated steady-state polarization pss (black solid lines) 
agrees reasonably with the experimental data^i. According to 
our analytical expression Eq. ®, the negative dip at B « 105 
mT and positive shoulder at B ~ 100 mT correspond, re¬ 
spectively, to IF- resonance at A = A/i//2 and IF+ resonance 
at A = -AnI2. The strong transverse HFI leads to rapid 
DNP and significant non-Markovian oscillation in the real¬ 
time evolution of the nuclear polarization p{t) [solid lines in 
Fig- Etb)], beyond the description of our rate equation theory. 
Going successively closer to GSLAC, the DNP time shortens 
[Fig. 13b)] and finally saturates at ~ 1.6 ps, lower bounded 
by the NV optical initialization time Tc ~ 1.1 ps. The small 
difference between the black and orange solid lines in Fig. 
|3a) confirms that the far off-resonant excited state HFI has a 
negligible effect. Indeed, the exact numerical solution to Eq. 
I® including the excited state HFI shows that over the whole 
range of magnetic field under consideration, the excited state 
HFI alone polarizes the nuclear spin towards pss ~ 0.5 with 
a DNP time ~ 30-40 ps (much slower than the ground state 
HFI induced DNP), consistent with the estimation based on 
the Fermi golden rule. 

Next we demonstrate that our analytical expression Eq. ® 
is accurate as long as the Markovian condition IF « 1/Tc is 
satisfied, e.g., for magnetic field away from GSLAC or for '^C 
nuclei with weaker transverse HFI. For this purpose, we man¬ 
ually scale down Ag^xx and Ag^yy by a factor rj to decrease the 
DNP rate, but keep the longitudinal component A^ j.^ invariant. 
In Fig. [3la), the analytical nuclear polarization pss agrees well 
with the exact solution to Eq. Q for 7/ > 5. In Fig. [Ifb), when 
approaching VF± resonances, the analytical DNP rate from Eq. 

sharply peaks and significantly exceeds 1 /tc ~ (1.1 ps) 
while the exact solution saturates to ~ (1.6 ps)~^. Away from 
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FIG. 3. (a) Steady-state polarization of a first-shell nu¬ 

cleus: exact numerical solution to Eq. (H (solid lines) vs. analyt¬ 
ical formula Eq. (dotted line). To illustrate the convergence of 
our theory, we scale down the transverse HFI (A^ j-j, and Ag^yy) hy 
a factor p = \ (black line), 5 (red line), and 25 (gray line), with 
Ag,xz - Ag,„ = 0 for T] = 5 and 25 to avoid spurious behaviors, (b) 
DNP rate for p = 25: exact solution to Eq. © (gray solid line) vs. 
analytical formula Eq. ([6]l (blue dotted line). 


W± resonances, our analytical DNP rate ^ 1 /tc and hence 
agrees well with the exact solution. 

In the above, for a first-shell nuclear spin, the separation 
- l-y^B - ~ -131 MHz is dominated by the strong 

HFI term. Consequently, the magnetic fields for the two res¬ 
onances at A = ±AnI2 are separated by A^HyeB) » 5 mT, 
which determines the magnetic held sensitivity of the polar¬ 
ization sign reversal in Fig. Oa). This suggests that for weakly 
coupled nuclear spins with much smaller « -2 MHz, 
the magnetic held sensitivity can be greatly enhanced. 


DNP of weak coupled nucleus 


Now based on the analytical expression Eq. ©, we pro¬ 
vide a complete picture for the dependence of the DNP rate W 
and steady-state polarization pss on the resonance linewidth 
r ^ R, the |0g)-| - Ig) separation A = Dgs - jgB, and the hy- 
perhne tensor Ag. It reveals the possibility to strongly polarize 
a weakly coupled, distant nucleus. 

The simplest case is strong optical pumping with the res¬ 
onance linewidth F > |A;vL i-e., the W± resonance peaks at 
A = -AnI2 and A = Av/2, respectively, are not spectrally re¬ 
solved: W± ^ n{\Ag^+zf.'^IA)6^^HA). In this case pss is uniquely 
determined by the HFI tensor Ag, but W still depends on 
the |0g)-| - Ig) separation A and becomes suppressed when 
I A| » F. For dipolar HFI, 


^g,ZZ — — — ^dp(|R|)(l + ^ 0 )? 

Ag,++=Adp(|R|)(l-ce)e2'^ 


where Adp(|R|) = -/io>riv/(47r|R|3) ^ 20 MHz AVlRp, 7n = 
-10.705 MHz/T is the nuclear gyromagnetic ratio, cg = 
3 cos^ 0-2, and 9 (if) is the polar (azimuth) angle of the 
location R. Thus the polarization 


Cg + 1 leg 


(7) 


is uniquely determined by the polar angle 9 and is independent 





FIG. 4. Steady state nuclear polarization pss [(a) and (b)] and DNP 
rate W in units of Hz [(c) and (d)] by dipolar HFI with the NV center, 
calculated from the analytical formula Eq. nil with optical pumping 
rates R = 3.2 MHz [(a) and (c)] and 0.2 MHz [(b) and (d)]. The 
distance of the nucleus from the NV center is |R| = 12 A in (c) 
and |R| = 20 A in (d). The white (black) dotted line indicates the 
critical magnetic field for WA (IV_) resonance. 


of |R| as long as the DNP rate W » 2ydep or equivalently 

A^pdRI) » 4Fydep, (8) 

e.g., forF = 0.2 MHz and ydep = 1 s \ this condition amounts 
to Adp > 1 kHz and is satished by nuclei with |R| < 27 A. An 
example for strong optical pumping is shown in Figs. |4ja) 
and life). In Fig. life), the W+ resonance peak (marked by 
white dotted lines) and the W- resonance peak (marked by 
black dotted lines) are signihcantly broadened along the mag¬ 
netic field axis and are not clearly resolved. Consequently, the 
nuclear polarization in Fig. Ufa) follows the HFI coefficients 
|Ag,+Tp oc (1 + cg)^ [plotted on top of Figs. Ufa) and|4fb)] and 
is well described by Eq. Q. 

The most interesting case is weak optical pumping with 
the resonance linewidth F ^ |AAf|. In this case, the W± res¬ 
onance peaks are spectrally resolved and can be selectively 
driven into resonance to achieve strong nuclear polarization 
Pss = =*=A^/(A^ -I- 2F^) ai +100% by choosing A = +A-[^I2. 
The DNP rate W cx 1 /F is dominated by the resonant process 
and can be increased up to the saturation value ~ 1 /tc by de¬ 
creasing the linewidth F. An example is given in Figs. |4fb) 
andllfd). In Fig. Hfd), the very narrow W+ resonance peak 
at A = -An 12 (corresponding to B « 102.37 mT, marked by 
white dotted lines) and the W- resonance peak at A « An/2 
(corresponding to B « 102.45 mT, marked by black dotted 
lines) can be clearly resolved. In Fig. Hfb), the W± res¬ 
onances give rise to fine structures in the nuclear polariza¬ 
tion pss near B » 102.37 mT and B 102.45 mT, super¬ 
imposed on the usual dependence on the polar angle 9 via 
W± oc |Ag,+:fP OC (1 + cg)^ [plotted on top of Figs. Ufa) and 
|4fb)]. We notice that a nucleus located at |R| = 20 A 
away from the NV center can be highly polarized with a DNP 
rate ~ 60 Hz = (2.6 ms) sufficient to overcome the slow 
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nuclear depolarization ydep ~ 1 s“'. 

Finally we present the three conditions for achieving strong 
nuclear polarization; 

1. The linewidth F « IAat] = \2yNB - Ag^^\, where ly^B ^ 
-2.2 MHz since we always consider B ~ 102 mT near the 
GSLAC. Therefore, IAat] ranges from ~ 2 MHz (dominated 
by ly^B) for weakly coupled nuclei to ~ 100 MHz (domi¬ 
nated by Agj^) for strongly coupled nuclei. 

2. For optimal performance, the |-lg)-|0g) separation A should 
be tuned accurately to match one resonance peak. This re¬ 
quires the mismatch 6A < F, |Aa^|, i.e., the magnetic field 
control precision 6B should be smaller than lA^^I/ye (~ 0.8 
G for weakly coupled nuclei with |Aa?| ~ 2.2 MHz near 
GSLAC) and F/y, (~ 0.07 G for F = 0.2 MHz), accessible 
by cutTent experimental techniques, e.g., 6B = 0.02 G has 
been reported^. If the control precision 6A » linewdith F, 
then the DNP rate would be reduced by a factor (dA/F)^. 
Taking |4|d) as an example, a large magnetic field control 
error SB = 1 G and hence dA = 2.8 MHz would reduce the 
maximal DNP rate from ~ 60 Hz to ~ 2 s \ which is still 
comparable with the depolarization rate ydep ~ 1 s '. 

3. For DNP to dominate over nuclear depolarization, the trans¬ 
verse HFI strength should satisfy 

-H » 4Fydep, 

which reduces to Eq. ([8]) for dipolar HFI. 


DNP of multiple nuclei 


Having established a completely picture for single-'^C 
DNP, now we generalize the above results to many (A » 1) 
nuclear spins coupled to the NV center via Hap = F, ■ 
I,, where F, = Sg ■ A, and the excited state HFI is neglected 
because the induced NV-nuclei flip-flop processes are off- 
resonant. In the same spirit as that in treating the DNP of a sin¬ 
gle nuclear spin, we decompose the HFI into the longitudinal 
Fi,zii,z the transverse part Zi^i(A,+7i,- + h.c.)/2. 
Then we approximate the longitudinal HFI with 
(which is treated exactly) and treat the transverse HFI with 
second-order perturbation theory, where 


^z — 'y ^ ^i,zz^i,z 


1=1 


is the collective Overhauser held from all the nuclei. 

The physics picture of the many-nuclei DNP is as follows. 
Up to leading order, the hip of different nuclei by the trans¬ 
verse HFI is independent, in the sense that at a given mo¬ 
ment, only one nuclear spin is being hipped, while other nu¬ 
clear spins simply act as “spectators”. However, the hip of 
each individual nuclear spin does depend on the states of 
all the nuclear spins via the collective Overhauser held hy. 
each many-nuclei state |m) = (|m,) is the Zeeman 

eigenstate of the /th nucleus) produces an Overhauser held 


hm = (mihjim) that shifts the energy of the NV state | - 1^) 
by an amount -h^- This renormalizes the | - lg)-|0g) sepa¬ 
ration from A = Dgs - yeB to A - and hence changes the 
NV dynamics and NV-induced nuclear spin hip, e.g., the spin 
hip rate lT,,±(hm) of the ith nucleus conditioned on the many- 
nuclei state being |m) is obtained from Eq. (|6ll by replacing 
Ag, A with A;, A - hm, respectively. 

The above physics picture is quantihed by the following 
rate equations for the diagonal elements = (in|p|m) of the 
many-nuclei density matrix p: 

Pm — ~ y [^i,+{^m)Pm ~ lF/^-(/Z|ji+ly)pm+l,] 
i 

- ^ [WiAhm)Pm - lT,-,+ (/jm-l,)Pm-lJ , (9) 


which has been derived by adiabatically decoupling the fast 
electron dynamics from the slow DNP process^^. Here 
|m -t-1;) = //"Im)/ is the same as |m) except that the state 
of the /th nucleus changes from |m,) to |m, -i- 1). Now we dis¬ 
cuss the difference between single-spin DNP and many-spin 
DNP. In the latter case, the DNP of each individual nucleus 
occurs in the presence of many “spectator” nuclei, which pro¬ 
duce a huctuating Overhauser held /r^ that randomly shifts the 
NV energy levels, such that the | - lg)-|0g) separation changes 
from A to a random value A - h^. This makes it more diffi¬ 
cult to tune the magnetic field to match the W+ (or VT_) res¬ 
onance exactly. More precisely, a finite mismatch ~ 
makes the originally resonant W+ (1T_) process off-resonant, 
and hence reduce the resonant DNP rate by a factor 
Eor example, a natural abundance of nuclei gives a typical 
Overhauser field {h^)rms ~ 0.2 MHz. This reduces the typical 
DNP rate by a factor of 2 for the linewidth F « = 0.2 MHz, 

e.g., the maximal DNP rate ~ 60 Hz [Fig. lUd)] for a 
nucleus at 2 nm away from the NV center is reduced to 30 
Hz » (5 ms) *, which is still sufficient to overcome the slow 
nuclear depolarization 7dep ~ 1 s“'. Therefore, under weak 
pumping F <si |A^|, the typical Overhauser field fluctuation 
(/Jz)rms ^ lA^I does not significantly influence the steady state 
nuclear polarization. 

For N nuclear spin-l/2’s, the number of variables (pm) is 
2^. When N is small, we can solve Eq. (|9]) exactly. Eor 
large N, however, the exponentially growing complexity pro¬ 
hibit an exact solution to Eq. (|9]l. Here, we introduce an 
approximate solution. The essential idea is to assume small 
inter-spin cori'elation, so that the A-nuclei density matrix p 
can be factorized as p - ®^_^pi, where the /th nuclear spin 
state Pi ^ I T/XT,' 1(1 + p!)/2 + \ |(1 - p,)/2 is de¬ 

scribed by its polarization p,. This dramatically reduces the 
number of variables from 2^ many-nuclei populations (pm) to 
N single-spin polarizations [p,). Substituting this approxima¬ 
tion into Eq. (|9]) and tracing over all the nuclei except for the 
/th nucleus give N coupled equations: 


Pi = 


-(W-+ + W,--) 



- W,.- \ 

+ W,-J 


( 10 ) 


where VT,',± is the spin flip rate of the /th nucleus averaged over 
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Magnetic field (mT) Magnetic field (mT) 



FIG. 5. Average polarization of (a) A = 7 and (b) N = 400 
randomly chosen nuclei for different optical pumping rates R = 
0.2, 0.8, and 3.2 MHz. For R = 0.2 MHz, the spatial distrihution 
of the polarization of A = 400 randomly chosen nuclei in a 
magnetic field (c) 102.36 mT and (d) 102.45 mT [marked by the 
arrows in panel (b)]. The depolarization rate ydep = 1 Hz for all 
panels. 


the states of all the other (N - 1) nuclear spins: 

Wi,+ = TrVF,-+(/ij)| i/Xi/ I pj, 

= TrWiXhz)\ T/XT/ I Pj, 


clearly demonstrate the possibility to strongly polarize weakly 
coupled nuclear spins by using weak optical pumping near 
W+ resonance (B » 102.36 mT) or W- resonance {B » 102.45 
mT). The steady state polarization of these weakly coupled 
nuclei is ultimated limited by nuclear spin depolarization. 
This can be clearly seen in the spatial distribution of the nu¬ 
clear polarization at the W+ resonance [Fig. |4|c)] and the W- 
resonance [Fig.|4|d)]. At W+ resonance, strong positive polar¬ 
ization is achieved for nuclei with |R| < 15 A away from 
the NV center, corresponding to dipolar HFI strength > 6 kHz. 
For more distant nuclei, the HFI strength is too weak for 
the DNP to dominate over the nuclear depolarization, so their 
polarization drops significantly. Similarly, at W- resonance, 
strong negative polarization can be achieved for nuclei 
with |R| < 25 A away from the NV center, corresponding to 
dipolar HFI strength ~ 1 kHz. In principle, further decrease of 
the optical pumping rate R and hence the resonance linewidth 
F could make the W± resonances even sharper, which allows 
more distant nuclei to be strongly polarized. 

Finally, we recall that the on-site or nucleus has a 
isotropic transverse HFI and hence W- - |Ag_++|^ = 0, i.e., 
they could be completely polarized and thus does not signifi¬ 
cantly influence the DNP of the nuclei, except for a shift 
of the I - lg)-|0g) separation A by ~ 2.2 MHz (for '"^N) or 
~ 3 MHz (for ^^N), corresponding to a shift of the resonance 
magnetic field by ~ 0.78 G (for ^^N) or ~ 1.1 G (for '^N). 


DISCUSSION 


i.e., Wi^± depends on the polarizations {pj} {j + i) of all 
the other nuclear spins. Therefore, Eq. (ITOl i with i - 
1,2, ■■■ ,N form N coupled differential equations for {/?,). 
The steady-state solutions {/7,,ss) are obtained by solving N 
coupled nonlinear equations with recursive methods. This 
approach provides a good approximation to the average po¬ 
larization pss - (1/A^) 2^1 Pi.ss and average Overhauser 
field /jss = but does not include any spin- 

spin correlation effect, such as feedback induced spin bath 
narrowingi^i^^^— . 

As shown in Fig. [Sja), for a small number (A = 7) of 
randomly chosen nuclei coupled to the NV via dipolar 
HFI, the average polarization pss from the above approxima¬ 
tion agrees well with the exact solution to the rate equations. 
For N - 400 randomly chosen nuclei, the exact solution 
is no longer available and we plot the approximation results 
in Fig. |3b). In both Figs. |5la) and 13b), we can clearly 
see that increasing the optical pumping strength significantly 
decreases the average nuclear polarization and its magnetic 
field dependence due to the broadening of W± resonances. For 
weak optical pumping R - 0.2 MHz, the average polarization 
Pss of 400 nuclei shows a positive maximum —1-25% at 
B » 102.36 mT (due to W+ resonance) and a negative maxi¬ 
mum -55% at Z? » 102.45 mT (due to W- resonance). The 

stronger polarization along -z axis arises from the dependence 
of the single-spin flip rate W± oc oc (1 + cg)^ on the po¬ 

lar angle 0 of the nuclear spin location: more nuclear spins 
have |Ag_++p > |Aj_+_|^ [as plotted on top of Figs. |4|a) and 
mb)] and hence favors negative polarization. These results 


In conclusion, we have presented a comprehensive theo¬ 
retical understanding on the dynamic nuclear polarization in¬ 
duced by an optically pumped NV center near the ground state 
anticrossing at ambient temperature. Our results not only pro¬ 
vide a clearly physics picture for a recently observed^ mag¬ 
netic field dependence of the polarization of first-shell nu¬ 
clei, but also reveals an efficient scheme to strongly polarize 
weakly coupled nuclear spins ~ 25 A away from the NV 
center (HFI strength ~ 1 kHz) by weak optical pumping and 
fine-tuning the magnetic field. These results provide a clear 
guidance for optimizing future dynamic nuclear polarization 
experiments. For example, this scheme could be used to polar¬ 
ize distant nuclei in isotope purified diamond^ to further 
prolong the NV coherence time. 


METHODS 

Numerical evaluation of nuclear spin transition rates 

Taking Wm+i^m as an example, we need to first calculate 
Pm from X.m,mPm = 0 and TrNvFm = 1, and then calcu¬ 
late {-Lm+\,m - ijNBy^F-Pm- For this purpose, we map NV 
operators into vectors and Liouville superoperators into ma¬ 
trices by introducing the complete basis set (l/y) s |/){ 7 |), 
where i,j - 1,2, • ■ • , M, and M = 7 is the number of NV 
energy levels in our model [see Fig. [Ha)]- Then the NV 
operator = Z/jUjlv) mapped to a vector v with 
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components V ,7 = {i\F-Pm\j), and the Liouville superoper- Here is the energy difference between |y)|m+1) and |/)|m) 
ator -Cm+i,m “ iJnS is mapped to a matrix S via - (|i), ly) are NV states and \m},\m ± 1) are nuclear Zeeman 

iyNB)\ij) - Y,ki \kl)Ski,ij, where the matrix element states), i.e., 


Ski.ij = {k\[(£m+i,m - iyNB)\i){j\]\l). 

Then we immediately obtain (-£m+i,m - iyNBy^F-P„, - 
Yjij |iy)(S“*v ),7 and hence the transition rates Wm+i<^m- 


Explicit analytical expressions for nuclear spin transition rates 


With I + Ig) and | + !«) neglected, the rotating frame Hamil¬ 
tonian of the five-level NV model [Fig. [Hb)] is 


//nv = ^O‘-lj,-lg+(£>es-7'cfi+W0“t^)o'-l„-l, + — 

Neglecting F^-induced NV spin mixing, the longitudinal HFI 
reduces to -H Ae^zzd'~i,-iPIz- Neglecting the 

non-collinear term oc Szl± (as the NV mostly stays in |0g)) 
and the transverse NV excited state HFI (which is far off- 
resonant), the transverse HFI reduces to (P+I- + h.c.)l2 
with F+ = (d--i^,o,Ag,++ -H d-Q^ _i^Ag _+)/ V 2 . To calculate 
from Eq. (01), we first determine the NV steady 
state as P,„^^ = ^g^o„o,, + (1 - ■Pg)o-o.,o,, where R = 
- to) is the optical pumping rate for 
|0g) <-> |0c). Substituting into Eq. (01) gives 

( 11 ) 

where is the (i, j) matrix element of the opera¬ 
tor p**^ s {£m±i,m + which is a linear combina¬ 
tion of and since the superop¬ 

erator £m±i,m for the five-level NV model [Eig. [Hb)] has no 
coherent coupling between {| 0 g), | 0 e)) and {| - Ig), | - l^)). 

Now we calculatep*^j^ g by taking the {0g| • | - lg),(-le| • 
| 0 g),{-lg| • | 0 e) and {-le| • | 0 e) matrix elements of 

(X„,±i,m + yNB)p^'^^ = 
which gives four coupled equations 


+ -Y^P-u, 




(A' 


(±) 

u,o. 




Q 




\ - I 

-U, 

A±) 


(A' 


(±) 

i„o, 


- 'r-i„o.)pit.o 


Qb 


Vto -P-to) = 0 - 


A_i^ 0,, ~ TeB ± y^B (jti ± l)Ag^^^, 

A i^ Op — f^es yeB i): yMB (ktl + l)Agzzr’ 

+ "0 - w. 

A-/^,o, = 'E-t.o, - "0 + W, 


and Fji is the linewidth of the NV transition |/) |y), i.e., 

T-ij.Oj = 7 ( 0 , r_4,o, = 7-1-71/2, F_i^,o, ^{Fe + y + yi +y^)l2, 

and F_i^,o, ^ {Te + y + y^)l2. Eliminating andp^*/^_g_^ 

gives the “rate equations”; 


from which we obtain the solution 


P-h.% 


^(±) 

-u.o. 


-/F . 0 -??W(l + 




A**’ 

-le.Oe 


-iT- 




( 12 ) 

where with = (QR/ 2 ) 2 /(Ai*\g^ - 

*F-i^,o.) and = {Q.Rl2f - iF.^^Oj) the complex 

self-energy corrections to |0g) and | - Ig) by the optical pump¬ 
ing |0g) «-> |0e) and | - Ig) I - 1?), respectively. More 
precisely, the optical pumping |0g) «-> lOg) (I - Ig) I - Ic)) 
induces an optical Stark shift Rey?[(^^ (-Re??*!*^) and an ef¬ 
fective dissipation Im'Rg'^^ (Im??^*^) for the NV ground state 
|0g) (I “ Ig))- Taking as an example, if |F_i^ oj is much 
smaller than g |, then the optical Stark shift 6Eq^ « 
{Q.rI2)^ !reduces to the conventional form of a second- 
order energy correction in non-degenerate perturbation theroy, 
while the dissipation takes the semiclassical form of a 

Fermi golden rule. Substituting Eq. (IT^ into Eq. (ITTl ) im¬ 
mediately gives VTm±i<_m, which assumes a tedious form as it 
includes various quantum coherence effects. 

Eor simplification, we use the fact that the NV excited state 
dephasing F,, ~ 10“* GHz » other NV dissipation y, 71 , 7 ,^ and 
typical detuning |wo - to\, g |, and g |, and further 
restrict to weak optical pumping R <si F_i^_o, (~ 26.3 MHz). 
In this case, the optical Stark shift is negligible and the optical 
pumping rate simplifies to R =s Q^/F^, so the self-energies 
=s iR only induces NV level broadening. Substituting the 
resulting expression g =s g - 1 ( 7 ^ - 1 - R)) into Eq. 

(ITT]) gives Eq. (|5]). 
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